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THE CORE-CUSP PROBLEM IN COLD DARK MATTER HALOS AND SUPERNOVA FEEDBACK: EFFECTS 
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ABSTRACT 

We investigate the dynamical response of dark matter halos against recurrent starbursts in forming 
less-massive galaxies to solve the core-cusp problem, which is a discrepancy between the observation 
and the cold dark matter model. The gas heated by supernova feedbacks after a starburst expands, 
and then the star formation terminates. This expanding gas loses energy by radiative cooling, and 
then falls back toward the galactic center. Subsequently, a starburst arises again. This cycle of 
expansion and contraction of the interstellar gas leads to the recursive change in the gravitational 
potential of the interstellar gas. The resonance between dark matter particles and the density wave 
excited by the oscillating potential plays a key role to understand the physical mechanism of the 
cusp-core transition of dark matter halos. The dark matter halos effectively gain the kinetic energy 
from the energy transfer driven by the resonance between particles and the density waves. We find the 
critical condition for the cusp-core transition that the oscillation period of the gas potential should be 
approximately the same as the local dynamical time of the dark matter halo. We present the resultant 
core radius of the dark matter halo after the cusp-core transition induced by the resonance using the 
conventional mass-density profile predicted by the cold dark matter models. Moreover, we verified 
the analytical model using TV-body simulations and the results nicely confirm the resonance model. 
Subject headings: galaxies: dwarf — galaxies: halos — galaxies: kinematics and dynamics — galaxies: 
structure — galaxies: evolution 



1. INTRODUCTION 

Cold dark matter (CDM) cosmology, which is the stan- 
dard paradigm of structure formation in the universe, 
contains several serious unsolved problems. Recent ob- 
servations of nearby less- massive galaxies and low surface 
brightness galaxies have revealed that the density profile 
of the dark matter halo is constant at the center of these 
galaxies (e.g., Moore 1994; Burkert 1995; de Blok et al. 
2001; Swaters et al. 2003; Weldrake et al. 2003; Gentile 
et al. 2004; Spekkens et al. 2005; Kuzio de Naray et al. 
2006, 2008; Oh et al. 2008, 2011). However, in cosmo- 
logical A-body simulations based on collisionless CDM, 
dark matter (DM) halos always have a steep power-law 
mass-density distribution at the center (e.g., Navarro et 
al. 1997; Fukushige & Makino 1997; Moore et al. 1998, 
1999; Jing & Suto 2000; Klypin et al. 2001; Fukushige 
et al. 2004; Navarro et al. 2004; Diemand et al. 2004; 
Reed et al. 2005; Stadel et al. 2009; Navarro et al. 2010; 
Ishiyama et al. 2011). This discrepancy is a well-known 
unsolved problem in the CDM scenario and is so-called 
the "Core-Cusp problem" . 

To solve this problem, this study has focused how 
baryon components gravitationally affect DM halos be- 
cause less-massive galaxies like dwarf galaxies are more 
sensitive to stellar activity such as supernova explosions 
than giant ones owing to their shallow gravitational po- 
tential well. In dwarf galaxies, the hydrodynamic re- 
sponse against the stellar activity depends on the combi- 
nation of various factors such as mass, morphology, and 
the star formation rate of the dwarf galaxies (e.g., Dekel 
& Silk 1986; Yoshii & Arimoto 1987; Mori et al. 1997, 
1999, 2002; MacLow & Ferrara 1999; Silich & Tenorio- 
Tagle 2001, Bland-Hawthorn et al. 2011). 
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In less massive galaxies with a smooth gas distribution, 
supernova feedbacks blow the interstellar gas out from 
the galaxy centers. This mass loss reduces the depth of 
the gravitational potential around the center of DM ha- 
los and may flatten the central cusp of these halos. So 
far, the dynamical effects of this process on DM halos 
have been studied using collisionless A-body simulations 
(Navarro et al. 1996a; Gnedin & Zhao 2002; Read & 
Gilmore 2005; Ogiya & Mori 2011; Ragone-Figueroa et 
al. 2012). Ogiya & Mori (2011) have argued the relation 
between dynamical response of DM halos and timescale 
of mass loss (i.e. the star formation rate). The den- 
sity profile of DM halos is correlated with the mass loss 
timescale and is flatter when mass loss occurs over a short 
time than that over a long time. However, even if mass 
loss occurs instantaneously, which is maximal effect, the 
central cusp remains; thus Ogiya & Mori (2011) conclude 
that the mass loss is not a prime mechanism to flatten 
the central cusp. When the mass loss timescale is much 
longer than the dynamical time of the system, the density 
profile of the DM halo recovers its initial cusp profile. 

Massive galaxies with moderate stellar activity are 
more robust against supernova feedback, so supernovae 
can expel only small fraction of the gas from these galax- 
ies. The gas heated by supernova feedbacks shortly af- 
ter a burst of star formation expands temporarily, and 
then the star formation is terminated. The expanding 
gas loses energy by radiative cooling and then falls back 
toward the galactic center. Then, star formation will 
enhance again and ignite star burst subsequently (Stin- 
son et al. 2007). This repetitive gas motion leads to 
the recursive change in the gravitational potential. Us- 
ing hydrodynamic simulations combined with star for- 
mation and supernova feedback, the dynamical effects 
of such recurring changes have been studied (Governato 
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et al. 2010; Pontzen & Governato 2012; Maccio' et al. 
2012) and indicate that oscillations in the gravitational 
potential may modify the central cusp into a flat core. 
Read & Gilmore (2005) also argue the importance of re- 
curring changes in the gravitational potential. The os- 
cillation timescale is also an important factor in deter- 
mining the dynamical effects on DM halos, and dwarf 
galaxies in the Local Group each have distinctive star 
formation histories (SFHs) (e.g., Mateo 1998; Tolstoy et 
al. 2009; McQuinn et al. 2010a, 2010b; Weisz et al. 
2011). Furthermore, recurring star formation activities 
duplicates properties (e.g., metal distribution, and the 
velocity-dispersion profile) of dwarf spheroidal galaxies 
(Ikuta & Arimoto 2002; Marcolini et al. 2006, 2008). 
On the other hands, several other mechanisms have also 
been proposed to solve the Core-Cusp problem. Repre- 
sentative models involve the dynamical friction induced 
by the motion of gas or stellar clumps (Mashchenko et 
al. 2006, 2008; Inoue & Saitoh 2011) and bar instability 
(Weinberg & Katz 2002). 

In this paper, we investigate the physical mecha- 
nism about the cusp-core transition under the recurring 
change in gravitational potential of the interstellar gas. 
Especially, we focus on the energy transfer driven by the 
resonance between DM particles and density waves in- 
duced by the oscillating gravitational potential. Then, 
we clarify the cusp-core transition of the DM halos using 
collisionless ./V-body simulations. This paper is struc- 
tured as follows. Our analytic resonance model is de- 
scribed in §2. In §3, we present the result of our numer- 
ical simulations and interpret them with the aid of the 
resonance model. Finally, we summarize and discuss the 
results in §4. 

2. ANALYTIC MODEL 

The oscillation process (i.e. dynamical oscillation of 
gas in galactic potential) is an external force-like forced 
oscillation for DM halos. What is the dynamical response 
of DM halos to this force? To answer this question, we 
construct a simple analytic model to analyze the pos- 
sibility of resonance between DM particles and density 
waves induced by the external baryon potential. 

2.1. Resonance model 

We approximate DM halo is treated as a fluid and the 
dynamical evolution of the system is described by the 
linearized fluid equations of the perturbed equilibrium 
state. Hereafter, we label physical quantities in the ini- 
tial equilibrium state as "0" , external forces as "ex" , and 
induced quantities as "ind" . We assume infinite constant 
density field to as the equilibrium state (pa = const.), 
pick out a group of particles which have some constant 
velocity (v — const.). In addition, the absolute values of 
the induced quantities are assumed to be much smaller 
than that of the equilibrium state or the external force. 
With these assumptions, the linearized equation of mo- 
tion is given by, 



dv ind (t,x) dv ind (t,x) _ d$ cx (t,x) 



dt 



v 



dx 



dx 



(1) 



where the induced potential $ lnd is neglected. 
Hereafter, we discuss the case in which the system is 



forced oscillation such as 



Asm (kx — fit), 



(2) 



where A, k and fl are the oscillation strength which is a 
positive constant, the wavenumber, and the angular fre- 
quency of the external force, respectively. The solution 
of Eq. (fTJ is given as follows, 



v ind (t,x) 



A 



fl — kvo 



cos (kx — fit). 



(3) 



The linearized equation of continuity can be similarly 
obtained, 



dp ind (t,x) ^ dp ind (t,x) 
dt dx 
and the solution is given by 



Po- 



dv i nd (t,x) 
dx 



Pind(t,x) 



Ap k 



(fl — kvo)' 



■ cos (kx — fit). 



(4) 



(5) 



From these solutions, we notice that when fl — > kvo, 
the coefficients diverge, and this is the resonant condi- 
tion. Unfortunately, the signs of induced values are not 
able to be decided without the relations between fl and 
kvo, and the phase of the wave. Next, we consider the 
energy interchange between the DM collisionless system 
and the oscillatory external force. 

In the following, we evaluate the change of the kinetic 
energy density of the particle group with the constant 
density p$ and velocity "Uo in the equilibrium state. The 
kinetic energy density in the equilibrium state Kq and 
after wobbled by the external force K are written by 

Ka = ^PoVo, (6) 

K = ~(po + Pind)(Wo + 2u Wind + vf nd ). (7) 

Thus the change of the kinetic energy density, SK is 

SK = K-K = ^po(2vov ind +vf nd )+^p iDd (vl+2voVi nd +vf nd ) 

(8) 

Averaging SK over a wavelength, we obtain the averaged 
change of the kinetic energy density, 



< SK >= 



p A 2 fl + kvo 
4 (fl - kv f 



(9) 



From Eq. ©, when fl — kvo > 0, the particle group gains 
the kinetic energy from the external force. On the other 
hand, when fl — kvQ < 0, its kinetic energy decreases. 

Note that we have thus far confined the argument to 
a particular group of particles. Practically, the system 
have uncounted groups of particles, and each of them 
has respective densities and velocities in the equilibrium 
state. In the realistic self-gravitational systems, particles 
having low energy (i.e. low velocity) are more dominated 
than particles having high velocity. Thus, the change of 
the kinetic energy density overall the system, AK will 
be positive, i.e., 



AK 



< SK > dp Q dvQ > 0. 



(10) 
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Fig. 1. — The integration path. The cross symbol shows the 
singular point, v = tt/k. We separate the contour into three parts, 
(i), (ii), and (iii). 



From Eq. (fTU)) . we know the system gains the kinetic 
energy from the external force substantively, and it will 
expand. In other words, v m d > and p ln d < for the 
whole of the system. 

As an example, we assume Maxwell-Boltzmann distri- 
bution for particle velocity and evaluate the change of ki- 
netic energy density overall the system, AK cx f °° dv < 

SK(v) > e~" 2 = J* °° dvf(v). To avoid the singular point 
of < SK(v) > (i.e., v = fl/k) we change the integration 
path as shown in FigJTJ The value of the integration is 
conserved through this change. The integration can be 
separated into three parts, 



AKoz 



(0 



dvf(v) 
dvf(v) 4 



dvf(v) 



(iii) 



dvf(v). (11) 



We can evaluate the second term using the residue theo- 
rem, and we find the second term is an imaginary num- 
ber. We ignore it. The first and the third are positive 
and negative, since Q — kv > and Q — kv < is satis- 
fied, respectively. The absolute value of the third term 
will be smaller than the first term, because it is damped 
by Gaussian strongly. Thus, AK will be positive, i.e. the 
system gains kinetic energy from the external force sub- 
stantively, and it will expand. In other words, Ui, lc i > 
and pind < for the whole of the system. 

According to the above argument, we conclude that 
the particles gain the kinetic energy effectively, if the 
particles satisfy the resonance condition, 



kv ~ fi. 



(12) 



Therefore, the system will expand to settle new equilib- 
rium state, and the density will decrease. 

2.2. Resonance condition 

In this section, we apply the resonance condition to a 
spherical system and simply extend it as follows, 



k(r)a(r) w ft, 



(13) 



where r is the distance from the center of the system, the 
wave number, k{r) as a function of r is defined by 



k(r) ee ^. 
r 



(14) 



Assuming the virial theorem, vq in Eq. (|12p is substi- 
tuted by the velocity dispersion 



cr(r) 



GM(r) 



(15) 



where M{r) is mass within r. Here, we define the dy- 
namical time of the system as 



where 



*d(r) = 



p(r) 



3?r 



Y 32Gp(r) ' 

M(r) 
(47r/3)r 3 ' 



(16) 



(17) 



is the average density interior to radius r. In this case, 
the resonance condition in Eq. (fT3| becomes 



*d(r) 



2^/2 



(18) 



where T ee 27r/ft is the oscillation period. 

Furthermore, we consider the resonances of the over- 
tone modes. Hereafter, "1" denotes the fundamental- 
tone, and "n" denotes n-th overtone (n = 2, 3, • • •). The 
resonance scale of the fundamental-tone is r\ = 2ir/ki, 
thus those of the overtones are 



2ir/k n = 2tt jnk\ — T\jn. 



(19) 



The resonance condition for the n-th overtone is td(r n ) ~ 
T/n. According to this argument, the resonance scale 
will be nearly equal to r\ for which the resonance condi- 
tion, 

*d(ri) ~ T, (20) 

is satisfied. From this equation, we can derive the reso- 
nant radius using the inversion procedure 



n = t- d \T) 



(21) 



The resonance of the overtone modes will occur in the 
central dense region of the DM halo and the resonance of 
the fundamental-tone will appear in the outer less dense 
region. In the next section, we apply this argument to 
the models of CDM halos predicted by cosmological N- 
body simulations, Eq. (|2"2")l and estimate the core scale. 

2.3. Core scale 

The mass-density profile of a DM halo obtained from 
cosmological TV-body simulations based on the CDM sce- 
nario can be fitted by the following formula: 



Pdm(x) 



Po 



x a (x ■ 



(22) 



where x = r/RoMi r is the distance from the center 
of the DM halo, a is a power-law index of the central 
cusp, and po and -Rdm are the scale density and scale 
length of a DM halo, respectively. This formula indicates 
that the density distribution of a DM halo changes from 
p oc r~ a in the center (r < -Rdm) to p cx r -3 at the 
outskirts (r > Rum)- Here, a = 1.0 corresponds to 
the Navarro-Frenk- White (NFW) model (Navarro et al. 
1996b; Navarro et al. 1997) and a = 1.5 is similar to the 
Fukushige-Makino-Moore (FMM) model (Fukushige & 
Makino 1997; Moore et al. 1999). 
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The mass profile derived by Eq. (|22[) is given by 

M(a; x) = 4ir f oJ?DM x 3 - a 2 F 1 [3 - a, 3 - a, 4 - a; -x] , 
3 — a 

(23) 

where 2 -F\ [3 — a, 3 — a, 4 — a; — x] is Gauss's hypergeo- 
metric function, which we denote hereafter simply by 2F1 
[a;-x\. From Eq.®, 



3 — a 



Po 



Km 



where the virial mass of a DM halo is defined by 



4"7T 

M vlr = yp cri (l + z) 3 Ai?: 



3 

vir ! 



(24) 



(25) 



where p cr j is the critical density of the universe, z is red- 
shift, A(= 200) is a parameter of the density enhance- 
ment, c = i? v i r /i?DM is a concentration parameter, and 
i? v i r is a virial radius. 

Using Eqs. (|2"3"|) and (|2~4")| . the average density of the 
DM halo interior to radius x is given by 



p(x) 



3M vir x~ a 2 Fi[a;-x] 



47ri?g, M c 3 - Q 2^1 [«; -c] 



(26) 



and by using the resonance condition given in Eq. (|20[) , 
we finally arrive at the core scale r corc = x corc -Rdm ■ 

We now discuss the utmost limit, x coro -C 1. When 
a^corc *C 1, the approximation 2*1 [a;— x] ~ 1 holds. 
Thus, the core scale is expressed as 



■ R 



DM 



8G 



ivJ vir ± 



1 !/« 



7T 2 R 3 



c 3 - a 2F![a; 



Rdm 



32G 
3tt 



Ap cri (l + z) 3 



(27) 



2*1 N -c] 



l/c 



(28) 



Using Eq. (|25[) . we obtain more convenient formulae 
for NFW model (a = 1.0) and FMM model (a = 1.5). 
The general formulae are given by 

1 ■> / M V i T \ 1/3 / T \2 

— e = 0.74^41; -a]" (1 + ,) (_) pc, (29) 

for the NFW model, and 

1/1 I M„ ir \ 1/3 / T \4/3 

r ra = 23 a *[1.6; + ,)(_^) (^) PC(30) 

for the FMM model. Here, we show some values of 
2-Fi[a; — c] in Tabled] for a = 1.0, 1.5 and some ordinary 
c parameters obtained from cosmological iV-body sim- 
ulations (e.g., Maccio' et al. 2008). When we assume 
c = 10 and z — 0, we reach the core scale by following 
simpler formulae: 



= 25 



M vlr a/3 



-10 9 M Q ; 
for the NFW model, and 

r -120 1 Mvir ^ 1/3 ' 

' corn — -Lavj 



10V 



pc, 



T 



,4/3 



pc, 



(31) 



(32) 



1O 9 M ; v 10 7 yr 

for the FMM model. 

In the next section, we compare our analytic prediction 
with the results of TV-body simulations. 



TABLE 1 

Numerical values of 2-Pi[a;— c] for conventional models of 

DM HALOS 



a (Model) 


c= 5 


c = 10 


c = 15 


c = 20 


1.0 (NFW) 


0.07667 


0.02978 


0.01631 


0.01046 


1.5 (FMM) 


0.16948 


0.08681 


0.05656 


0.04117 



Note. — a is the power-law index of the central cusp of a DM 
halo and c is the concentration parameter. 



3. NUMERICAL SIMULATIONS 
3.1. Numerical models 

In this section, we describe our collisionless TV-body 
simulations. To generate TV-body systems with a cusp 
described in Eq. (|22[) , in the equilibrium state, we use 
the fitting formulation of the distribution function (DF) 
proposed by Widrow (2000). In this case, the DF only 
depends on energy, thus the velocity dispersion of the 
system is isotropic. 

We adopt the Hernquist potential (Hernquist 1990) to 
describe the baryon potential around the center of the 
DM halos. The external potential is given by 



$ b (r,t) = — 



GMh 



R b {t) : 



(33) 



where Mb and Rb(t) are total baryon mass and scale 
length of the external potential, respectively. We change 
the scale length of the baryon potential to Rb(t) 



0.5[R K 



Rb,imn] [1 + co s (^)] + #b,min, where tt 



2tt/T,T is the timescale of oscillation of the external po- 



tential, Ri 



and Ry, 



are maximum and minimum 



scale length, respectively. 

To construct a stable initial condition, we dynamically 
relax the equilibrium TV-body system quoted above in 
the external baryon potential with the fixed scale length 
i?b = -Rb.max to ~ 30 id (0.2), where id (0.2) is the dy- 
namical time defined by Eq. (fT6"| at r = 0.2 kpc. The 
resultant density profiles are almost the same as the ini- 
tial non-perturbed system of NFW and FMM model. 

In this study, we simulate the dynamical response of 
a DM halo with virial mass M v j r = 10 9 M Q , the virial 
radius i? v h- = 10 kpc, the scale length -Rdm = 2 kpc, 
and we assume z — 1. In this case, the mean dynami- 
cal time within 0.2 kpc is t d (0.2) = 4 Myr for the FMM 
model and 10 Myr for the NFW model. Other param- 
eters adopted in each of our simulation runs are shown 
in Table [2j Throughout this paper, the total number of 
particles is N = 1,048,576, the tolerance parameter of 
the Barnes-Hut tree algorithm (Barnes & Hut 1986) is 
9 = 0.8 and the softening parameter is e = 0.008 kpc. 

3.2. Results 

In this section, we present the results of our iV-body 
simulations and compare them with the theoretical pre- 
diction via the simple analytic model given in §2. 

Figure [2] shows the distribution function of DM par- 
ticles as a function of the radial velocity. Dashed line 
is the equilibrium (non-perturbed) state of the iV-body 
system for FMM model. Before adding the oscillatory 
external potential, we generate an initial stable TV-body 
system with a fixed external baryon potential for a con- 
stant scale length Rb(t) = i?b,max defined by Eq. ([33f . 
Thin dotted line corresponds the distribution function of 
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TABLE 2 

Summary of simulation runs 



ID 


T[t d ] 


M b [ xlO s M Q ] 


-Rb.max [kpc] 


-Rb.min [ 


1'MMl 


i 


1.7 


2.0 


0.04 


FMM2 


3 


1.7 


2.0 


0.04 


FMM3 


10 


1.7 


2.0 


0.04 


FMM4 


0.1 


1.7 


2.0 


0.04 


FMM5 


50 


1.7 


2.0 


0.04 


FMM6 


3 


1.7 


1.0 


0.04 


FMM7 


3 


1.7 


2.0 


0.1 


FMM8 


3 


1.0 


2.0 


0.04 


NFW1 


1 


1.7 


2.0 


0.04 


NFW2 


3 


1.7 


2.0 


0.04 


NFW3 


10 


1.7 


2.0 


0.04 




0.1 1 o.i 

rlkpc] 



Fig. 3. — Density profile (top panels) and spectrum ol radial 
velocity (bottom panels) of the DM halos after several oscilla- 
tions of the external potential. Left panels are results of FMM 
model (FMM1, FMM2 and FMM3), and right panels are NFW 
model (NFW1, NFW2 and NFW3). Solid, dotted, dot-dashed, 
and dashed lines represent T = 1, 3, andlOtd and the initial condi- 
tion, respectively. Each vertical line is the core scale predicted by 
Eq. 03 and J5D). 

panels correspond the core scale predicted by Eq. (|29|) 
or ([30]) . where the left and right lines are the cases for 
T = i<j (0.2) and T = 3id(0.2), respectively. As you can 
see, our analytic model perfectly predicts the core scale 
derived by these numerical experiments. On the other 
hand, the system evolves adiabatically for the slow oscil- 
lation with T = 10td(0.2). In this case of the adiabatic 
change in potential, the system always accommodates 
the potential, and the central cusp remains (Ogiya & 
Mori 2011). Although the central cusp is remains, how- 
ever, a bump appears at the outskirts of the DM halo. 
This bump is obviously generated by the resonance of the 
density wave. Despite beyond the appropriate range of 
the approximation a; coro -C 1 , Eq. (|30|) gives a reasonable 
prediction such as r corc = 0.95 kpc for FMM halo. 

Here, we compute the Fourier transform of the radial 
velocity of the systems as a function of r using each snap- 
shot of data to verify the emergence of the energy transfer 
due to the resonance between the particles and the den- 
sity waves. In the lower panels of Fig. [3j we show the 
Fourier components of the radial velocity of the angular 
frequency, ui, equal to the oscillation frequency of the ex- 
ternal potential Q. As you can see, there are the peaks 
accord with the core scale in the density profile, and the 
position of these peaks depends on T. We have demon- 
strated the spectrum of the fundamental mode for each 
run. These results clearly show that, as expected, the 
core scale agrees with the position of the peak predicted 
by Eq. ([2TJ|) . As we discussed in the previous section, 
the resonance of the slow oscillation appears at outer 
low-dense regions of DM halo, whereas that of the rapid 
oscillation appears at the inner more dense regions. This 
result indicates that when DM halos are in resonance 
with the external potential, they effectively accelerate 
particles, so that the core is created around the position 
of the peak. 

Figure [4] shows the spatial distribution of the Fourier 
components of the radial velocity including the case of 
u> ^ CI. Although we only show the results of the FMM 
model in this figure, a similar phenomenon appears for 
the NFW model. In the bottom panels of Fig. [3J we 
only show the fundamental mode for each T. Figure [4] 
clears up the resonance of overtone modes. For exam- 



Note. — FMM and NFW indicate initial DM halo model. T, 
Mb, Rb max an d Rb min are the period of oscillation, the total 
mass, and the maximum and the minimum scale length of the 
external baryon potential, respectively. We adopt the mass ratio 
of baryon component to a DM halo cosmic value obtained from 
Wilkinson Microwave Anisotropy Probe observations (Spergel 
et al. 2007; Komatsu et al. 2009, 2011). 
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Fig. 2. — Distribution functions of DM particles for FMM model 
as a function of the radial velocity. Dashed line is the equilibrium 
(non-perturbed) state. Thin-dotted line (thick-dotted line) is for 
t = 15td(0.2) (t = 30td(0.2)) after adding the oscillatory exter- 
nal potential with a fixed external baryon pot enti al for a constant 
scale length Rb (t) = Rb max defined by Eq. (133 I t . Solid line cor- 
responds the distribution function for the quasi-equilibrium state 
after several oscillations. 

DM particles at 15t d (0.2). After 30£ d (0.2) , the system 
reached the quasi-stable state, and it remains virtually 
unchanged (thick dotted line). Then, we add the oscilla- 
tory external potential given by Eq. (|33[) . Solid line cor- 
responds the distribution function with the system after 
several oscillations. It clearly shows that the number of 
low- velocity particles decreases and that of high- velocity 
particles increases. This indicates that the system gains 
the net energy through the acceleration of the slow par- 
ticles by the density waves induced by the oscillatory 
baryon potential. 

Figure [3] demonstrates how the dynamical response of 
DM halos depends on the oscillation frequency of the 
external baryon potential. In each panel, solid, dotted, 
dot-dashed and dashed lines correspond to T = 1, 3 and 
10£d(0.2) and the initial condition, respectively. The 
upper-left (upper-right) panel corresponds the final den- 
sity profiles of DM halos for initially FMM (NFW) model 
after several oscillation periods. All of the models clearly 
show the cusp-core transition and the resultant core 
scale depends on the oscillation frequency of the exter- 
nal baryon potential. The vertical dashed lines in both 
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Fig. 5. — Density profile of the DM halos after oscillations of 
the external potential. Dashed line is the initial condition of the 
FMM model. Dotted line and thin-solid line represent T = O.ltd 
(FMM4) and T = 50t d (FMM5), respectively. The dot-dashed line 
expresses the maximum-contracted case Rb(t) = 0.04kpc using the 
adiabatic contraction discussed by Blumenthal et al. (1986). 



Fig. 4. — Spectrum of radial velocity of the DM halos (FMM1, 
FMM2 and FMM3). Each line represents the same for Fig. |3]and 
each panel corresponds for (a) T = t&, (b) T = 2t d , (c) T = 3t d , 
(d) T = 4t d , (e) T = 5i d , and (f) T = Wt d . Vertical lines show the 
predicted resonance scale, r n . We consider the fundamental mode 
to correspond to T = 10t d (viz. u>\ = 27r/10td). 

pie, the dot-dashed line, T — 10id(0.2), exhibits peaks 
in panels (f), (e), (b), and (a). These peaks correspond 
to the resonance of the fundamental mode and overtone 
modes n = 2,5 and 10, respectively. Furthermore, verti- 
cal lines indicate the resonance scale between the funda- 
mental mode and the overtone modes predicted by Eq. 
(Tl9|) . If we can measure a resonant scale, then we can 
predict other resonance scales using this relation. 

Figure [5] shows the resultant density profiles of the 
FMM model for the cases of the extremely rapid oscilla- 
tion (model of FMM4) and the extremely slow oscillation 
(model of FMM5). Dotted line and thin solid line repre- 
sent T = O.H d (0.2) and 50i d (0.2), respectively. Dashed 
line is the initial condition, and dot-dashed line expresses 
the maximum-contracted case with Rb(t) = 0.04 kpc us- 
ing the adiabatic contraction model argued in Blumen- 
thal et al. (1986). For the rapid oscillation model, the 
DM particles are not able to respond the rapid change of 
the external potential. Therefore, the energy exchange 
between the density waves and the DM particles is in- 
efficient, and the central cusp remains. In other words, 
this is limited by a numerical resolution. In the high- 
resolution simulations, the resonance will appear at the 
central region because the resonance is a scale-free. It 
should be noted that DM particles are exposed to some 
fixed external potential (see dotted line at the center), 
the DM halo concentrates more than the initial FMM 
model, and the degree of concentration is weaker than 
that for the maximum contraction case (dot-dashed line). 
The solid line corresponds the nearly adiabatic change in 
the gravitational potential. In this discussed by 

Ogiya & Mori (2011), the system always accommodates 
the changing potential, and the central cusp remains, al- 
though there seems to be a bump in the mass density in 
the outskirts. The scale of the bump is nearly five times 
that for T = I0t d (cf. Fig. E}. 

Figure [6] demonstrates how the dynamical response of 




r[kpc] 

Fig. 6. — Density profile of the DM halos after oscillations of the 
external potential. The thin solid line is the initial condition, FMM 
model. Thick solid, dotted, dot-dashed, and dashed lines represent 
the result of FMM2, FMM6, FMM7, and FMM8, respectively. 

DM halos depends on the other parameters (Mb,i?b,max 
and i?b,min). The thin solid line is the initial condition 
of the FMM model. Thick solid, dotted, dot-dashed and 
dashed lines represent the resultant density profile of the 
model of FMM2, FMM6, FMM7, and FMM8, respec- 
tively. It is apparent that the resultant core scale is 
almost alike for each model, and the final density pro- 
file of the simulation is quite similar. Our results show 
that these parameters affect the number of oscillations 
needed to reach quasi-equilibrium state with the central 
core, but there is almost no effect on the final state of 
DM halos. This indicates that the period of oscillation 
is the most important parameter to determine whether 
the central cusp is flattened, and the core scale. 

4. SUMMARY AND DISCUSSION 

We have studied the dynamical response of DM halos 
to recurring changes of the gravitational potential of the 
interstellar gas around the center of the DM halos. The 
resonance model between particles and the oscillating ex- 
ternal potential is proposed to understand the physical 
mechanism of the cusp-core transition of DM halos. It is 
found that the collisionless system effectively gains the 
kinetic energy from the energy transfer driven by the res- 
onance between DM particles and density waves induced 
by the oscillation of the gravitational potential of the 
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TABLE 3 

Comparison with observational results 



H 



NGC 2366 
NGC 6822 
rg II 



[kpc] p(rcore) [M R pc- J ] t d ( 



[My] 



[Myr] 



0. t [ 



0.05 
0.03 
0.03 



SO 



100 
100 
100 



Note. — Approxii 
core observed cor. 



lated values of NGC 2366, NGC 6822, and 
radius, p(T- corc ) is mass-density at T- cor c , re 
obtained these values from Wcldrakc et al. (2003; NGC 6822), 
(2008, 2011; NGC 2366, Holmberg II). We estimated the dynamical 



ilmbcrg n. 
itively. We 
Oh et al. 



*d( r co 
obta" 



e ) using Eq. 1161 . tq F is the into 
id from McQuinn et al. (2010a, 2010b). 



al of sta 



fori 



interstellar gas. The condition for the cusp-core transi- 
tion is that the oscillation period of the baryon potential 
should be approximately the same as the local dynami- 
cal time. In other words, the resonance condition at r is 
represented by id (f core) = T and the overtone modes are 
also in the resonant states. In addition, the core radius 
of the DM halo after the cusp-core transition driven by 
the resonance between particles and the density waves is 
shown using the conventional mass-density profile of DM 
halos, which is predicted by the cosmological structure 
formation models. Moreover, we verified our analytical 
model using TV-body simulations and the simulation re- 
sults nicely confirm our resonance model. Therefore, we 
conclude that the resonance between the DM particles 
and the oscillation of the baryon potential plays a key 
role in the cusp-core problem of DM halos. 

Pontzen & Governato (2012) have a similar approach 
to understand how recurring potential flatten out the 
cusp of DM halos. They have studied about two lim- 
ited cases in which potential changes instantaneously or 
adiabatically, and found that most particles in the sys- 
tem gain energy when potential change occurs instanta- 
neously, and systems expand. As seen from above, our 
resonance model gives broad understanding of the fun- 
damental physical mechanism of the cusp-core transition 



of DM halos. In addition, the model is able to predict 
the core scale quite precisely. 

The oscillation of gas is directly correlated with the 
stellar activity of galaxies. In recent years, SFH and 
the density profile of nearby dwarf galaxies have been 
studied precisely. Our results imply that these two phe- 
nomena are correlated. For example, we adduce three 
galaxies; NGC 2366, NGC 6822, and Holmberg II (Table 
[3]) , which have a core structure at the center ( Weldrake et 
al. 2003; Oh et al. 2008, 2011) and episodic SFHs (Mc- 
Quinn et al. 2010a, 2010b). Their dynamical time at 
the core scale corresponds approximately to the interval 
of their episodic star formation activity. These galaxies 
conceivably justify the resonance condition predicted by 
Eq. ([20]) at the core scale. Yoshii & Arimoto (1991) 
have studied the color change of galaxies with oscillating 
star formation history. They revealed that the colors of 
galaxies also change oscillatory and that star formation 
activity imprints its record, depending on the initial mass 
function and oscillation period, in the two-color diagram. 
Evidence of oscillating SFHs may be detected, and the 
oscillation period might be measured more precisely in 
the future observations. 

In the series of subsequent studies, we will examine a 
more realistic model that includes non-spherical analysis 
and gas dynamics and that will more closely study the 
correlation between density structure and the SFH of 
dwarf galaxies. 
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Grants-in-Aid for Scientific Research: (A) (21244013), 
(C) (18540242), and (S) (20224002), and Grants-in-Aid 
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